Theory of spin-polarized scanning tunneling microscopy applied to local spins 
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We provide a theory for scanning tunneling microscopy and spectroscopy using a spin-polarized 
tip. It it shown that the tunnehng conductance can be partitioned into three separate contributions, 
a background conductance which is independent of the local spin, a dynamical conductance which 
is proportional to the local spin moment, and a conductance which is proportional to the noise spec- 
trum of the local spin interactions. The presented theory is applicable to setups with magnetic tip 
and substrate in non-coUinear arrangement, as well as for non-magnetic situations. The partitioning 
of the tunneling current suggests a possibility to extract the total spin moment of the local spin 
from the dynamical conductance. The dynamical conductance suggests a possibility to generate 
very high frequency spin-dependent ac currents and/or voltages. We also propose a measurement 
of the dynamical conductance that can be used to determine the character of the effective exchange 
interaction between individual spins in clusters. The third contribution to the tunneling current is 
associated with the spin-spin correlations induced by the exchange interaction between the local spin 
moment and the tunneling electrons. We demonstrate how this term can be used in the analysis of 
spin excitations recorded in conductance measurements. Finally, we propose to use spin-polarized 
scanning tunneling microscopy for detailed studies of the spin excitation spectrum. 



PACS numbers: 72.25.-b,73.63.-b 

I. INTRODUCTION 

Pushing the limits for detection of electronic, mag- 
netic, and vibrational properties towards the quantum 
limit requires appropriate experimental tools and tech- 
niques. In particular, single atomic spinJ^^ and niag- 
netic nanostructures consisting of few magnetic atompnSl 
on nonmagnetic substrates are frequently studied as 
model systems for miniature data storage devices, spin- 
tronics apphcations, and qubits which are crucial for 
quantum information technology. Being well-defined and 
controllable on the atomic scale, they ideally serve their 
purpose for studies of fundamentals of their local prop- 
erties and interactions. 

The scanning tunneling microscope (STM) was in- 
vented in the 1980'sPfor the purpose of imaging metallic 
surfaces with atomic resolution, and it was theoretically 
described by TersofF and Hamann,'^^! relating the tunnel- 
ing conductance to the density of states (DOS) of the 
local environment. More than being a scanning tool, one 
can park the tip over an object and perform differential 
conductance measurements in order to reveal informa- 
tion of the local electronic structure. Using the STM 
equipment as a means for spectroscopical measurements 
was previously discussed, see e.g. Refs. fTTHTOl the ap- 
proach that is commonly known as scanning tunneling 
spectroscopy (STS). This technique has since been suc- 
cessfully applied in several directions, e.g. detection of 
noise,!— indirect measurements of the Kondo effect,!^ and 
the observation of exchange splitting.ESl 

Extensions of the STM/STS techniques have been pro- 
vided by using a spin-polarized tip (SP-STM/SP-STS), 
e.g. magnetic Cr02 tips,'^^' or e.g nonmagnetic W tips 



which are coated with a ferromagnetic metal, e.g. 
or an anti-ferromagnetic metal, e.g. To this end, 

the theory by TersofF and Hamann was extended to 
also account for the spin-polarization in the tip and 
substrate. f^SHsU The latter formulation was for instance 
used in Ref. "JB" in analyzing experimental SP-STM re- 
sults of spin-moments for single Co atoms on a Pt sub- 
strate. The lack of explicit reference to the spin moments 
of the adatoms located on the substrate in this theory, 
however, calls for an advancement in the theoretical for- 
mulation. In this paper we present a theory which sug- 
gest to use the SP-STM technique for directly extracting 
quantitative information about the magnetic moments of 
nano-scale objects. 

The theory presented here is discussed in the context 
of recent experimental observations, using both STM and 
STS approaches and using both a nonmagnetic and spin- 
polarized tip. Our discussion will be cast in the light 
of the theoretical description of the SP-STM, particu- 
larly for measurements performed in presence of local 
spin moments, S„, located at r„ on the substrate. We 
find that indeed STM tunneling is in principle capable 
of detecting single spin and moreover detect the spin ori- 
entation. This conclusion can be drawn, for instance, 
when we consider the tunneling electrons to be interact- 
ing with the local spins through exchange. For this type 
of interaction mechanism, the tunneling matrix element 
can be separated into one spin-independent and one spin- 
dependent component. Under such conditions, we find 
that the tunneling conductance can be separated into 
three components, of which the first provides a conduc- 
tance depending on the electron and magnetic densities 
of the tip and substrate, the second yields a conductance 
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which is directly proportional to the local spin moment, 
and the third being proportional to the noise produced 
by the local spin fluctuations. The last contribution to 
the current, or conductance, was recently discussed in 
Refs. MM 

We point out that we are interested in the qualitative 
effects caused by the presence of local spin, or magnetic, 
moments located on the substrate surface. For this rea- 
son, we will build upon previous results concerning the 
matrix elements for the tunneling electrons between the 
tip and substrate. In particular, the spin-dependence in 
the tunneling current generated by the local spin mo- 
ment via exchange interaction is of interest, while other 
features that are pertinent to the tunneling, e.g. geom- 
etry of the tip, surface, and adsorbate, bias voltage de- 
pendence of the matrix element, etc, will be treated with 
less accuracy. 

We begin the paper by a derivation of the tunneling 
current and the corresponding (differential) conductance 
in Sec. [n] We continue by discussing the main properties 
of the conductance contributions in Sec. Ill and there- 



after we discuss possible experimental measurements sug- 
gested by the theoretical results in Sec. |lVj and we con- 
clude the paper in Sec. [Vj 



II. PROBING THE LOCAL SPIN MOMENT 

The wave function of the tunneling electrons in the 
tip, separated from the substrate by distance d, has an 
exponentially small overlap with the substrate electron 
wave functions. The spin-dependent tunneling matrix 
element can be calculated using Bardeen's resultjl^ and 
is given by^i^ 



r„ =ro exp - 




$- JS„(i)-cr 



$0 



(1) 



which is understood as a matrix in spin space. Here, 
(T = (cr^, cr^, cr^) is the Pauli matrix vector, whereas Tq 
describes the spin-independent tunneling in absence of J. 
$ is the tunneling barrier height, while 'I'o = fi^/SmcP 
is the energy related to the distance d between the tip 
and surface. In principle, the spin-dependent tunnel- 
ing matrix element depends on energy, position, ap- 
plied bias voltage, and quantum numbers of the elec- 
tron stat es, which have been discussed extensively in the 
literaturcfiSElIllHlZl Jn the present study we shall omit 
those dependencies, however, for the sake of focusing on 
the dependence of the tunneling conductance on the local 
spin moments. 

The exchange energy J|S| is small compared to the 
barrier height, so that we can expand the exponent and 
find the effective tunneling matrix element 



where 

To.roe-V^coshfl yj:, (3a) 
T,.roe-V^sinh:f (3b) 

such that Ti/Tq ^ J/$. For metals and semiconductors 
it is reasonable to use J ^ 0.1 eV ,123 while the tunneling 
barrier $ ~ 1 eV, giving typical values of Ti/Tq ^ 1/10. 
In the following discussion, the tunneling rates Tq and Ti 
are treated as constants, except that we allow Ti to carry 
a spatial dependence which is related to the positions of 
the local spin moments. More details of this description 
of the tunneling matrix elements can be found in e.g. 
Refs. MM 

We next assume that the substrate surface is metallic 
for which "Hsub = J2ka ^ktr c^^Ckcr is sufficient, where cl^ 
creates a surface electron with energy Sko- , momentum k, 
and spin a. The energy-momentum dispersion relation 
need not be of free-electron character but may assume 
any general form, for which the specific details are unim- 
portant for the present derivation. The energies eko- are 
given relative to the Fermi level ep, which is common 
for the system as a whole. We associate the electronic 
and magnetic densities A^(r,e) and M(r, e), respectively, 
with the electrons in the substrate. For simplicity, we will 
assume in the following that those densities are slowly 
varying with energy. Analogously, we model the elec- 



trons in the tip by 



tip 



E 



per ^po-'^po-'-pcrj 



and define 



its corresponding electronic and magnetic densities n{e) 
and m(£), respectively. 

In general, the magnetic moments of the substrate and 
tip may be in a non-coUinear arrangement. Thus, defin- 
ing the z-direction of the global reference frame in e.g. 
the spin-quantization axis of the tip, the operators of the 
substrate are transformed according to 

Ckt \ cos(6l/2)e-*'^/2 -sin(6l/2)e-*'*/2 \ / Ck+ 
Cki y sin(6'/2)e^'^/2 cos(6i/2)e^'^/2 ) \ Ck- 

(4) 

where and 6 are the azimuthal and polar angles between 
the local substrate reference frame (spins s = ±) and the 
global one (spins a =t,i)- 

Tunneling of electrons between the tip and substrate 
in the presence of the local spin moments is captured by 
the model 

'^T=Y. ''lAToSaa' + Ti(r - r^O'T,.' • S„] 



pkn 



^Tq + Tier ■ Sr, 



(2) 



X Cko-'G 



where the bias voltages applied across the junction is 
denoted by eV. Due to the local nature of the spins, 
we have included a spatial dependence in the inter- 
acting tunneling rate, and we use e.g. Ti{r — r„) — 
T' exp (— |r — r„|/A), where A is the decay length. 



3 



The charge current running between the tip and the 
substrate is derived using non-equilibrium Green func- 
tions on the Keldysh contour, starting from the funda- 
mental relation I{t) = —edt^^^{n-p„) , where n^a = 



•!paCpa, givmg 



X e 



pkn 

era' 



(6) 



where we have defined T^^r' i) ~ ToScrcr' + Ti{r — 
T^n)o'(ja' ■ S„(i). The current expresses the rate of 
change of the expectation value of npo-, correspond- 
ing to the number of charges, e, in the tip. Expand- 
ing the average in the above expression according to 
{A{t)) « {-i) ]^{[A{t),'HT{t')])dt' , and by converting to 
real times, the current can be written as 



/(r,t,y)=|Re^^^/_* 



)-r+ieV (t~t') 



pp (Tcr 
kk' (t'V" 



cL'" it')f,'"<y" {rm,t')cp,,„ {t')])dt'. (7) 

Here, notice the appearance of the commutator taken 
between the operators inside the average. Separating 
the tunneling operator T into its components Tq and 
Ti, we find that the current can be naturally written 
as a sum of three terms, i.e. I{t) = J2i=o^''i^)- 
consider the differential conductance dI{r,t,V)/dV = 
Si=o 9Ii{r,t,V)/dV of stationary source drain voltages 
since such conditions are predominant in experimental 
situations. 



A. Background conductance 

The first contribution from the average in Eq. ([t]) is 
given by 

T^S,,, 5,.,.,. { [cl, (t)ck. (t) , 4,^ (i')cp'. it')]) 
=ro''^.."5."'."((4,Wcpv(t'))(ck.(i)4,.(<')> 



(cpv(t')4.W)(cL.(i')ck.W», 



(8) 



where we have decoupled the averages in the integral 
in terms of correlation functions for the substrate and 
tip electrons. Under the assumption that scattering 
between different momentum and spin states can be 
neglected, and introducing the notations g^„{t,t') — 

i{cLit')c^At)) and g>At,t') = {-i){c^At)cLit')), the 
first contribution to the tunneling current is defined by 



Ioir,t,V)=- 



2eT^ 



-Rc 



EE 

pk CT 



(9) 



Here, the lesser and greater correlation functions of the 
tip electrons are given by 



</>(t,i')=(±*)/(±£p.)e-'^-(*-*'\ (10) 



5, 



respectively, and analogously for the correlation func- 
tions of the substrate electrons. Here, f{x) is the Fermi 
function. Thus, performing the time-integration and re- 
placing the momentum summations by integrations over 
the corresponding spin-resolved density of electron states 
(DOS) ncr{e) and Na{r,e) in the tip and substrate, re- 
spectively, we finally arrive at the tunneling current 



Io{r,t,V) 



^ J n,{s-eV)N4r,e) 



X [/(£ - eV) - f{e)]de 
--—^ j[f{e)-f{e~eV)][n{e-eV) 
X N{v,e) + xa{e - eV) ■m{r,e))d£. (11) 



The last expression is obtained by noting that = (n + 
am^)/2 and A^^ = {N+aM^cos0)/2 (the factor a ±1) 
which thus provides a description that is independent of 
the coordinate system. 

We are, ultimately, interested in the tunneling (dif- 
ferential) conductance, dI{r,V)/dV, to which Iq con- 
tributes 



dV AksT J 2kBT 

X iV(r, e) -I- m(e - eV) ■ M(r, e)]de 

^TT^aoTg[n{eF - eV)N{r,eF) 

+ ni{sF - eV) ■M{r,eF)), T ^ 0, (12) 



where ctq = 2e^ /h is the fundamental conductance unit, 
whereas fcs and T is the Boltzmann constant and temper- 
ature, respectively. The expression is obtained under the 
assumption that the electronic and magnetic densities in 
the tip vary slowly with energy. While this conductance 
contribution relates the DOS and magnetic density in the 
tip {n{eF — eV), m{eF — eV)) and substrate {N{r,eF), 
M(r,ei?)), it is independent of the local spin moments. 



B. Dynamical conductance 

We next consider the expression for the first order con- 
ductance in the total local spin moment. This is provided 
by a similar consideration of the second contribution to 
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the tunneling current, i.e. beginning from 
Ji(r,t,y)=?^Re5]<,5]ri(r-r„) f [f{s^,) 

x/(-ek.)(S:(t')>-/(-£p.)/(£k.) 

X (5^(t))]e'(^>"'-^''<'+'^^)(*-*')dt'. (13) 

Here, we again have assumed that scattering between 
different momentum and spin channels within the tip 
and substrate is negligible. We replace the time- 
dependent spin by its Fourier transform, i.e. {S^[t)) ~ 
J (S',^(a;))e*"*dw/(27r). Along with the replacement of the 
momentum summations by energy integration over the 
corresponding spin resolved DOS of the tip and substrate, 
we find that this current can be written 



h{r,t,V) 



^Ti(r-r„) 1(5^ H) cos c.^ 



X [n{e)M^{Y,e')cose + m^{e)N{v,e')] 

X [f{e)f{~e')5{e-e' -io + eV) 

- f{-e)f{e')6{e - e' + eV)]dede''^. (14) 



We, thus, obtain the conductance 



dh{y,t,V) _7rVoro 



^Ti(r-r„) j{S^^{Lo))cos{ujt) 



dV AkeT 

X [m^{s)N{r,e') + n{e)M^{r,e') cos{0)] 

X [6{e -e' -uj + eV)+ S{e - e' + eV)] 

X cosh"^ — - — dsde''^. 

2fcsT 27r 

X N{r,eF) +n{e)M^{r,eF)cos{e) 
X [5{e + eV - uj) + 5{e + eV)]de^, 

ZTT 

(15) 



as T ^ 0. 

C. Conductance with spin-spin correlations 



The last term contained in the original expression given 
in Eq. ([7| can be written as 



/2(r,t,y) =^Re^^ri(r-r„)ri(r-r„) f n,{e) 

aa' nm 

J —OO 

X ■ (S„(i)S„,(t')) • T,'. - f{-e)f{e') 
X rT„,„ ■ (S„(i')S„(i)) • (T„„,]dt'dede\ 

(16) 



explicitly expressed in terms of the spin-spin correlation 
functions of the local spins. Although the spin-spin cor- 
relation function e.g. (S„(i)Sm(t')) provides all infor- 
mation from the local spin correlations to the tunneling 
current, it is convenient to rewrite this general function in 
terms of the propagators xi^^(t, i') = {—i){Sn{t)S^{i')) 
andx^m(t,i') = {~i){Sl{t)S^{t')). We notice in the sta- 
tionary regime, that these correlation functions depends 
on the time-difference t — t' , which allow us to express 
them through the Fourier transforms, e.g. X^m(^i^') ~ 
/ X^m(w)e-'"(*-*')da;/(27r). These remarks lead to that 
we can write the third contribution to the tunneling cur- 
rent according to 
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h{r,t, V) =^ ^Ti(r - r„)Ti(r - r„,) / [fie)fi-e')Sie - e' - u: + eV) - J{-e)J{e')5{e - e' + u + eV)] 



([X,t™(^) + XnU^)]He)N{e') - m(e) • M(e')] 



+ AxlMHe)N{e') + m{e) ■ M{e')] dede' 



m,{e)N{e')] 
duj 



2tt' 



(17) 



Despite the unappealing length of this expression, it provides a convenient starting point for analyses of spin inelastic 
tunneling spectroscopy. The conductance corresponding to this tunneling current becomedP^ 

^^^1^ (^) ' ^ E ^1 - - ^'») / (/(^)'^(^ -s'-uj + eV)+ f{-e)5{e - e' + c + eY) 

X {[xtM + XnMMs)N{s') - m(£) . M(e')] 
+ [X,t™(^) - X-+(w)][n(£)M,(£')cos0 - m,(e)7V(e')] 



+ 4x:^(w)[n(ef )7V(e) + m{e) ■ M{e')] cosh 



— ^ — dede' — 
2kBT 2tt 



^ ri(r - r„)Ti(r - r„) f (/(e)(5(e - c.- + eV) + f{^e)d{e + + eV) 

nm 

X ([x;tm(^) + Xni{^)]W{e)N{eF) - m(£) • M(e;^)] 
+ [Xn™(^) - x;;+(c^)][n(£)M,(£;.)cos0 - m, (e)iV(£j.)] 



+ 4x:„(w)[n(e)iV(£j^) + m(£) • M(£f )] )de^, T ^ 0. 



(18) 



We summarize this section by pointing out that Eqs. adatoms. The basic formulas, Eqs. (12), (15), and (18 1, 



(12 1, (151, and (181 serve as the main results of this pa- 



per. Below, we shall clarify the use of the partitioning of 
the conductance and how we can interpret experimental 
results in terms of the different contributions. 

We notice, however, that the electronic structures in 
the tip and substrate are described in a simple fashion, in 
that we assume dependence on a single momentum vector 
in terms of the homogeneous GFs, e.g. gka{t,t') for the 
substrate electrons. Beginning from the current given in 
Eq. ^ we can, however, repeat the above derivation al- 
lowing for an inhomogeneous description of the electronic 
structures, using e.g. .go-o-' (k, k'; i, t') for the substrate 
electrons. This will modify the formulas given for the 
conductance, something which goes beyond the scope of 
the present study. 

In connection to this, we also notice that we, here, have 
disregarded that the electronic and magnetic densities in 
the substrate are influenced by the local spin moments of 
the adatoms. In reality, the local spins and the substrate 
mutually influence one another and self-consistently cre- 
ate a total effective magnetic moment locally around the 



remain unchanged, however, by a self-consistent treat- 
ment of the substrate and adatoms. From the point of 
view we take in this paper, we can assume that the self- 
consistently created magnetic moment in the substrate is 
already included in the density M. 



III. DISCUSSION OF THE TUNNELING 
CONDUCTANCE 

The flrst contribution to the conductance, see Eq^ 



essentially captures the Tersoff-Hamann theorjESI^^ 
for magnetic tip and substrate. One should notice, how- 
ever, that we here have neglected the speciflc relation 
between the local tunneling matrix element Tq and the 
local DOS in the substrate, since we here are interested 
in the possibility to resolve its dependence on local spin 
moments located on the substrate. In studies of mag- 
netic surfaces, this conductance provides a sufficient tool 
for analysis of the magnetic structure of the surface, by 
using the spin-polarization of the electronic structure at 
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the Fermi level, since the density M(ei?) contains all such 
information. An analysis of local spin moments adsorbed 
onto the substrate can, in this model, only be performed 
on the level of the DOS and spin splitting of electron 
states at an energy defined by the Fermi level and the ap- 
plied bias voltage. The total spin moment of the adatom 
can, nevertheless, not be accessed through this conduc- 
tance. This fact is based on that although one can include 
the adsorbate into the electronic and magnetic densities 
of the surface, one only obtains energy dependent infor- 
mation of the adsorbate, while the spin moment is defined 
by the spin-polarization integrated over all energies. Any 
deeper information about the adsorbate can, thus, not be 
achieved. In this sense, we can consider the contribution 
dlo/dV as generating a back-ground conductance which 
includes features originating from the substrate. 



The second conductance contribution, see Eq. (15), 



contains terms which depend on the spin-polarization in 
the tip, nizie), and of the substrate, Mz{r,eF), respec- 
tively, and vanishes in case both the tip and substrate 
are non-magnetic. An important result here is that this 
contribution explicitly identifies a linear relationship be- 
tween local spin moment, (S'^(cj)) = J {S^{t))e~^'^*dt, 
and the tunneling conductance. Eq. ( 15 ) implies that 



one, in principle, can obtain a direct estimation of the size 
of the local spin moment simply by measuring the differ- 
ential conductance at its location. In order to achieve 
this functionality, it is necessary, as seen in Eq. (151, 



that the spin-polarizations of the tip and substrate are 
known. We further discuss this application below. It 
is worth pointing out, however, that although we have 
assumed stationary, i.e. time-independent, conditions, 
there is nevertheless a time-dependent component in the 
tunneling conductance. This time-dependence is gener- 
ated by the dynamics of the local spin and it is this fea- 
ture that would enable a read-out of the local spin mo- 
ment. The occurrence of the time-dependent component 
in the conductance, and tunneling current, also suggests 
a mechanism which may be employed for generation of 
high frequency electrical currents. The effective magnetic 
fields acting on the local spins, e.g. anisotropy fields etc., 
may be of the order of meV, see below, which would en- 
able generation of high GHz to THz ac currents. 



The last conductance contribution, see Eq. (18 1, pro- 



vides signatures that are generated by the spin-spin cor- 
relations, or spin fluctuations, occurring in the adsorbate. 
It is important to note that this conductance is finite for 
any polarization of the tip and the substrate, even when 
both electrodes are non-magnetic. Hence, regardless of 
the electronic and magnetic conditions of the system, this 
conductance directly depends on spin fluctuations in the 
adsorbate. It is also noticeable that the formulation of 
dl2/dV suggests one, in principle, will be able to study 
and distinguish between particular spin excitations in the 
adsorbate. This functionality is expressed from the fact 
that the sum and difference of the correlation functions 
Xnm Xnm multiplied by different combinations 
of the electronic and magnetic structures of the tip and 



substrate. One is therefore capable to configure the STM 
set-up in order to probe particular spin excitations in the 
local spin. This will be discussed further below. 



IV. PHYSICAL INFORMATION CONTAINED 
IN THE TUNNELING CONDUCTANCE 

We now discuss a few different physical examples which 
are introduced in order to shine some light on different 
aspects of the tunneling conductance. Before we go into 
the examples, however, we introduce the model of the 
local spins adsorbed on the surface. 

We consider a cluster of spins on the substrate and 
write the Hamiltonian for the spin S„ in the cluster 
according to = gfJ-B^ • S„, where B is an exter- 
nal magnetic field. The effective exchange interaction 
between the spin moments in the cluster is given by a 
Heisenberg model Hj — —JJ2n^m^n ■ S^. This ef- 
fective exchange comprise a combination of e.g. direct 
Heisenberg exchange and RKKY-like (Ruderman-Kittel- 
Kasuya-Yosida) exchange. The sign of the effective J 
may, thus, vary with distance between the spins in the 
cluster In this way we describe clusters of spins with a 
total spin moment S. We define the resulting eigensystem 
{Ect, \S,cr)}, a = —S, —S + 1, . . . ,S, for the eigenenergies 
and eigenstates, respectively, of the model 

n n 

(19) 

where we have added the spin anisotropy fields D and E. 
This model is pertinent to the recent studies of transi- 
tion metal elements adsorbed onto surfaces e.g. Fe and 
Mn on CuN surfacem,f^ and Fc and Co on Ft surface.'^ 
In those elements, the magnetism is constituted by the 
electrons in the d-shell whereas the electrons in the s- 
and p-shell strongly hybridize with the surface states in 
the substrate. We, therefore, assume that the s- and 
p-electrons in the adsorbate are included in the Hamilto- 
nian for the surface electrons and, below, we shall refer 
to them as conduction electrons. The electrons in the 
d-shell can be regarded as only weakly hybridizing with 
the conduction electrons which, therefore, enables us to 
disregard their contribution to the conductance. 

Under the assumption of very weak coupling between 
the d-electrons and the conduction electrons, which is 
pertine nt to the recent experiments on e.g. Fe and Mn 
on CuN,!^ we write the total dynamical local spin mo- 
ment in terms of the eigensystem defined for the cluster 
according to 



s 



(20) 



where f{£<j) provides the occupation of the state 
\S,a). Hence, the conductance is p rop ortional to 
S{S'{uj))dLo/{2n) = Y.^af{e^), c.f. Eq. (pll below. 
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In order to also describe the experiments conducted 
on e.g. a Pt surface^ we phenomenologically introduce 
a width of the local spin excitations by replacing Dirac 
delta functions by Lorentzian functions. Such a treat- 
ment can be physically motivated by that this would be 
the essential effect of a mean field approximated dressing 
of the bare spin averages and spin-spin correlation func- 
tions. This phenomenological model is sufficient for our 
present purposes since we do not attempt to explain all 
the details of the experiments. 



A. Estimating the spin moment 



The conductance in Eq. 



( 15 1 can be directly linked to 
This is simplest seen 



the total local spin moment (S^ 
by assuming that the electronic and magnetic densities in 
the tip and substrate vary slowly with the energy. Calcu- 
lating the Fourier transform of dli (t) jdV and integrating 
over all frequencies gives 



9/i(r,y) 
dV 



n 

X N{r,eF) +n{eF - eV)M:,{r,eF)cos{9) 

(21) 




25 50 75 
time (ps) 



FIG. 1: (Color online) Various aspects of dli/dV for a sin- 
gle 5" = 1 adatom. (a) dh(t,V)/dV, Eq. for B = Bz 
and mtip anti-parallel, B = — 1 T (lower plots), and par- 
allel, i? = 1 T (upper plots). The plots correspond to 
T = 0.3 K, = (dashed), T = 4.3 K, = (solid), 
(b) and (c) \dIi{u),V, B) / dV\ as function of uj (horizontal 
axis) and B (vertical axis) at T = 0.3 K, and (b) E = Q 
and (c) _E = 2 meV. Bright (dark) colors correspond to large 
(small) amplitude. The conductances have been normalized 
by 27rVoTori[m^Af + nM-,\. Here, D = -10 meV, » = 2, 
= 3n(£F)/4, = 0, and V = 0.3 meV 



where we have identified the total spin moment of the 
nth spin by {S^) = J {S^{uj))duj/{2n). In this fashion, 
we obtain a linear relationship between the differential 
conductance dl/dV and the average spin moment (S^). 

Our theory suggests that the STM conductance gener- 
ates a time-dependent signal, in agreement with earlier 
studies which suggest time-dependent noise spectroscopy 
using STM tunneling current, '^^ regardless of the time- 
dependence, or time-independence, of the bias voltage. 
The Fourier transform of this signal provides an energy 
resolved signal from which additional understanding, i.e. 
dynamics, about the spin systems can be extracted. 

Consider the simple example, given for e.g. a single 
(n = 1) spin moment of 5 = 1, for which a = 0, ±1. 
The anisotropy fields D and E break up the spin sym- 
metry, such that the spin levels are given by eg = and 
£±1 = D± ^/W+ljjfj^jBy , for B = Bz. The dynamical 
spin moment, thus, becomes {S^{uj)) — 7r[f(ei)S{uj—ei) — 
f{s-i)S(uj — £-i)], which generates a non-vanishing cur- 
rent when the populations /(e±i) of the states |1, ±1) are 
different. This observation holds for an adatom with any 
value of the spin and emphasizes the fact that the local 
spin has to have a definite moment in order to be mea- 
surable through the SP-STM. Using this dynamical spin 
moment in Eq. (15) and assuming slowly varying elec- 



tronic and magnetic densities in the tip and substrate, 
results for low temperatures in 



dhjt) 
dV 



■^/(si) cosEii — /(e_i) cose_ii. (22) 



More generally, the time-dependence of the conduc- 
tance can be written 



dhjt) 
dV 



a f{s„) cos Sat. 



(23) 



The period of the conductance oscillations are, thus, 
directly linked to the energy levels of the spin states, see 
Fig. [T] (a) for an example of the time-dependent conduc- 
tance for two different set-ups of the SP-STM. The period 
of the conductance oscillations can, thus, be changed by 
applying an external magnetic field Bz in order to vary 
the energy levels e^r , as is illustrated in Fig. [T] (b) , show- 
ing the magnetic field dependence of the Fourier trans- 
formed dynamical conductance. The time scale associ- 
ated with the conductance oscillations is given by the set 
of eigenenergies Eg., which means that only very low ly- 
ing excitation energies (0.066 meV corresponding to 100 
GHz) are reachable by means of the state-of-the-art ex- 
perimental technology. Although detection of spin mo- 
ments using this method definitely challenges todays ex- 
perimental resources and capabilities, it should be an ac- 
cessible regimes within the nearest future. 

It is important to notice, however, that the presence of 
the time-dependent component in the conductance opens 
the possibility to generate high frequency ac currents. 
Due to the high anisotropy fields acting on the local 
spin moment, of the order to meV (see Sec. IV B), one 



would be able to generate electrical and spin-dependent 
ac currents and/or ac voltages with frequencies in the 
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THz regime. 



C. Character of exchange interaction parameter 



B. Anisotropy parameters 

Our next example of the usefulness of Eq. ( |2l| ) is to 
consider an 5* = 1 adatom, and we plot for this sys- 
tem the conductance dli/dV, see Fig. [l](a), where the 
anisotropy fields D — —10 meV and E = meV for dif- 
ferent temperatures, i.e. T=0.3 K (dashed) and T=4.3 
K (solid). These para meter s are pertinent to the recent 
studies of Co/Pt(lll).'iiS126|The plots in (a) clearly show 
the phase shift of tt when the magnetic field B is reversed 
from anti-parallel (lower) to parallel (upper) orientation 
with rritip. The magnetization curves are S-shaped (not 
shown) , which correspond to a paramagnetic behavior of 
the local spin and which was observed in experiments.^^ 

From the dynamical conductance, dI{t)/dV or 
dI{Lu)/dV , one can obtain information about the mag- 
netic anisotropy fields D and E acting on the local spin. 
In the case of D < and E = 0, the spin moment points 
perpendicular to the surface, which generates a maxi- 
mal conductance. The response remains strong for val- 
ues of B sufficiently large to maintain a spin-splitting of 
the adatom which is larger than the thermal excitation 
energy ksT. This is illustrated in Fig. [l] (b), where 
\dil(u!, B) /dV\ is plotted for an = 1 adatom, as func- 
tion of u! and B. Bright (dark) colors correspond to large 
(small) amplitude of \dil{uj, B)/dV\. The lines that rep- 
resent the ground state of the adatom cross at zero mag- 
netic field, which is expected since the ground state is 
spin-degenerate at i? = 0. 

A finite value of £', on the other hand, provides a tilting 
of the spin moment away from the perpendicular orien- 
tation, which leads to a weakened conductance. This is 
illustrated in Fig. [l] (c), where E = 2. The finiteness 
of E breaks the spin-degeneracy of the ground state at 
vanishing magnetic field. The ground state is, however, 
pointing parallel to the substrate, which leads to that 
(S'^(cj)) = at i? = 0. The resulting magnetization line 
remains S-shaped (not shown). Despite the low temper- 
ature, however, the magnetization curve is reminiscent of 
the magnetization curve for higher temperature. This is 
expected since the anisotropy field E introduces an en- 
ergy barrier for the adatom to overcome in order to point 
its spin perpendicular to the substrate. 

In the production of the plots in Fig. [l] (b) and (c) we 
replaced the Dirac delta function by a Lorentzian func- 
tion with a broadening of 0.01 meV. This broadening is 
roughly 2 order of magnitudes smaller than what is ex- 
pected from the experimental set-up with Co/Pt(lll), 
see Sec. |IVD| for estimates, and is chosen in order to 
emphasize the different behavior of the conductance of 
different anisotropy parameters. Using more realistic 
anisotropy parameters blurs the resulting image, how- 
ever, the main difference between the character of the 
conductance for vanishing and finite field E can be re- 
solved. 



Studies of the dynamical conductance can also be used 
to reveal the sign of the effective exchange interaction 
between magnetic adatoms located on the surface. Con- 
sider for example a spin dimer. In case of ferromag- 
netic exchange, J > 0, the ground state of the spin 
dimer is a spin triplet. Varying the magnetic field, the 
ground state of the dimer acquires a magnetic moment 
e.g. {S^) = S or —S, depending on whether the mag- 
netic field is anti-parallel or parallel with the magnetic 
moment in the tip, in analogy with the single spin case 
discussed above. Hence, there is a measurable dynamical 
conductance. In case of an anti-ferromagnetic exchange, 
however, the ground state of the spin dimer is a spin sin- 
glet, with zero magnetic moment. Then, the dynamical 
conductance vanishes, c.f. Eq. (15), and the total con- 



ductance is strictly time- independent, except for possible 
noise fiuctuations. 



D. Spin fluctuations 

To illustrate the effect of spin fiuctuactions, we con- 
sider a local spin moment S comprising two coupled spins 
S„, n = 1, 2, a spin dimer, and consider them to be anti- 
ferromagnetically coupled. The ground state is a spin 
singlet \S = 0,(7 — 0), while the first excited states con- 
stitute a spin triplet jS" = + 82,(7 = 0, ±(5"! + 82))- 
Assuming an exchange energy | J| = \Es — Et\ > ksT, 
where i?T(S) denotes the triplet (singlet) energy, in order 
to prevent thermal excitations at zero bias, the equilib- 
rium conductance is given by the elastic tunneling be- 
tween the tip and the substrate only, i.e. dl/dV — 
dlo/dV. Effects from tunneling electrons scattering off 
the local spin moment averages to zero. 

The coupling to the tunneling electrons via the spin- 
spin interaction e.g. c|,^cro.cr' • S„Cko- enables, on the 
other hand, each individual spin constituting S to un- 
dergo spin-flip transitions which are assisted by spin-flips 
of the tunneling electrons. Due to this coupling, the cor- 
relation function e.g. o-g-cr' • (S„(i)Sm(t')) • cTo-'v is non- 
vanishing, in general. The spin-spin interaction, thus, 
provides a coupling between the singlet and triplet states 
which supports transitions between them. As a result of 
these transitions, a new channel for conductance opens 
at bias voltages V > \J\/e. 

The spin-spin correlation functions x^+(w), xtmi^)i 
and Xnmi'-^) calculated in terms of the eigensystem 
for the total spin, giving 



Anm 



{u)=i-^)2nJ2{^\8^\v){v\8^\^)PiE,) 

iv 

x[l-P{E,)]6iuj + E,-E,), (24a) 
x'nM=i-^)2nY.{^\S'Jv){v\8:J^)P{E,) 

iv 

x[l-P{E,)]6iuj + E,-E,), (24b) 
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FIG. 2: (Color online) dl/dV for Fe/CuN under different external magnetic fields at T = 0.5 K; (a) B = Bz, (b) B = Bx, and 
(c) for different F. Plots are off-set for clarity. 



where i (v) denotes the initial (intermediate) state, 
whereas P{x) accounts for the population in the corre- 
sponding state. Eq. (24 1 provides the general qualitative 



features of the spin-spin correlation function, and shows 
that there will appear steps in dl/dV whenever the bias 
voltage matches the transition energy zL{Ei — Ey). In the 
case of a spin dimer, for instance, there appear steps in 
dl/dV when the bias voltage supports the inelastic tran- 
sitions between the single and triplet states. In case of 
e.g. Fe/CuN, the calculated results are plotted in Fig. 2] 
for a non-magnetic tip and substrate, showing that the 
presented theory reproduces the results discussed in Ref. 
UHland shows an excellent agreement with experiments.!^ 
Signatures of the excitations in the experimental mea- 
surements do have a finite width, which corresponds to 
that the intermediate states have finite lifetimes. Replac- 



ing the delta functions in Eq. ( 24 ) by Lorentzian funtions 
l/(a:;^-|-F^), we phenomenologically include the (uniform) 
lifetime h/T for all intermediate states. In case of a single 
Fe on CuN with spin S = 2,^ and the anisotropy param- 
eters given in Tab. |Tj we find that the Fe spin is weakly 
coupled (F ~ 10 - 30 ^eV) to the Cu(lOO) through the 
CuN layer, see[2](c). This value is extracted by compar- 
ing the ratio between the maximal and minimal conduc- 
tance with the experimental result (~ 1/2). Similarly, 
we also extract the widths of the single Fe {S = 3/2) 
and Co {S = 1) adsorbed onto Pt(lll) surface.^ In Fig. 
[s] (a) and (b), we have plotted d'^I/dV'^ for Fe (upper) 
and Co (lower) and we find the best correspondence with 



TABLE I: Anisotropy parameters, D and E, used for the con- 
ductance plots given in Fig. [2] (c) andjs] the best widths, F, 
for the intermediate states in the spin-spin correlation func- 
tions. 



s 


D (moV) E (meV) 


r (meV) 


Fe/CuN 2 


-1.55 0.31 


0.01-0.03 


Fe/Pt(lll) 3/2 


-3.25 


5 


Co/Pt(lll) 1 


-10 2 
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experiments using the parameters given in [T] In Fig. [3] 
(c), wc finally provide a computation of the Fe dimer on 
Pt(lll) reported in Ref. iSj Because of the presence of 
two Fe atoms, we make use of space dependence of the 
tunneling rate and the decay length A, and plot the com- 
puted d^I/dV^ for different values of the decay length. 
In comparison with the experimental results, the decay 
length of 0.5 A gives the best agreement, suggesting a 
very rapid spatial decay of the spin-dependent tunneling 
rate. In the figure, we for reference have also included 
the resulting d^I/dV^ in absence of the broadening. 

In the calculations we have estimated the popula- 
tion numbers P{Eii^y)) of the states \i{v)) involved in 
the spin-spin correlation functions, c.f Eq. ( |24| ), by 
making use of the following observation. By expand- 
ing the correlation function e.g. Xram(iji') in terms 
of its eigenstates, we can write it as xtmi^^t') = 

where 4 

(di) creates (annihilates) a particle in the state \i). 
In the atomic limit, we can employ the decou- 
phng {{dld,){t){dldi){t')) = {dlm{t')}{d,{t)dl{t')} = 
P{E,)[l - P(£;^,)]e'(^*"-^-)(*"*'), where the population 
number P{Ei) = {d\di) can be estimated by using e.g. 
the Gibbs distribution P{E) = e-^^/Y.^ e"'^^'- 



E. Probing specific spin excitations 



(|18|) and (|24|) that the transi- 
v){v\S+\i}, and 



We observe in Eqs 
tion sequences {i\S^\v){vlS^\i), (Iffi 
{i\S^\v){v\Sjj^\i) are associated with different projections 
of the spin-resolved LDOS in the tip and substrate. 
The sum and difference of the first two sequences cou- 
ple to n{e)N{eF) — m(e) • M(ei?) and mz{s)N{eF) — 
n{e)Mz{ep) cos9, respectively, while the last sequence 
couples to n(e)iV(ei?) -l-m(e) •M(eF). The different cou- 
plings reflect an ability to enhance or attenuate the re- 
sponse of certain inelastic transitions, at will, by using 
different combinations of electronic and magnetic densi- 
ties in the tip and substrate. 
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In STM without spin-polarization (m, M = 0), for in- 
stance, the excitation spectrum can be analyzed to cer- 
tain detail by means of applying an external magnetic 
field e.g. along the z-direction, of the spin quantization 
axis of the sample. Such application of the magnetic 
field introduces a Zeeman splitting of the levels, which 
thus leads to a separation of the peaks in the (Pl/dV^. 
To be specific, consider a single adatom with 5 = 1, 
and anisotropy fields D < and E = 0, which is de- 
scribed by the states = l,cr = 0,±1), where E„ < Eq 
at B = 0. In this system, the leading contribution to 
the term in dl2/dV containing Xnmi^) + xtmi^)j c.f Eq. 
m, is proportional to E.=±i[^(^^-) " P{EQ)][f{eV ~ 
.iff + Eq) — f{eV + Ea - Eq)], while the term contain- 
ing Xnm{^) ~ Xnmi'-^) vanishcs. The former contribu- 
tion generates steps in the conductance at bias volt- 
ages eV = ±(£'cr — Ea). Those steps are separated by 
\E+i — whenever the states |1, 0) and |1, a) are un- 
equally occupied. Hence, the spin splitting imposed by 
the external magnetic field is detectable through sepa- 
rate but equally high steps in the conductance. This 
splitting can be seen in Fig. [4] where we plot d^I/dV^ of 
this scenario, where the lowermost curve correspond to 
non-magnetic conditions, whereas the second curve from 
below is obtained with B = Bz and B — 1. 

Using SP-STM opens further possibilities in the studies 
of spin systems, since then the term in dl2/dV containing 
— xtmi'^) is finite. In case of the S = 1 adatom, 
the leading contribution to this term is proportional 
to Ea^ii^lPiEa) + P{Eo) - 2P{E^)PiEo)][f{eV - 
E„ + Eq) — f{eV + E„ — Eq)], where it is important 
to notice that the contributions to this term have op- 
posite signs. Hence, by combining the electronic and 
magnetic densities in the tip and substrate such that 
TUzN — nM.zC0s6 < 0, this term leads to an attenu- 
ated (intensified) signal from the transition |1, 0) (1, -|-1| 




1 
bias voltage (mV) 



FIG. 4: (Color online) dPl/dV'^ for a spin S = 1 impurity, for 
various conditions of the external magnetic field B, magne- 
tizations of the tip (m) and substrate (M), and the relative 
angle (9) between the magnetization directions. Plots off-set 
for clarity. 



(|1, 0) (1, — 1|), c.f. third and fourth curves from below 
m Fig. [i) By the same token, the signal is intensified 
(attenuated) when m^A^ — rtM^ cos 6* > 0, see uppermost 
curve in Fig. |4] While the details of the terms containing 
the sum and difference of the spin-spin correlation func- 
tions Xnm ^'^d Xnm differ from system to system, the 
general conclusion we can draw out of this observation 
is that the intensity of the signal from any specific tran- 
sition depends on the magnetic densities of the tip and 
substrate, and on their relative orientation. 



F. Maximizing the signal from fluctuations 
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FIG. 3: (Color online) d^I/dV'^ for (a) Fe/Pt(lll) and (b) 
Co/Pt(lll) for different widths F, and (c) 2Fe/Pt(lll) for 
different decay lengths A, at T = 4.3 K. In panel (c), d^I/dV^ 
with zero broadening and A = 0.5 A is plotted for reference. 



Finally, we notice that the conductance dIo/dV van- 
ishes for nN -I- m • M = 0, which corresponds to the 
case with a half-metallic tip and substrate such that their 
magnetic moments are in anti-parallel alignment. In this 
set-up also the conductance dli/dV = 0, which implies 
that the measured signal is generated solely by dl2/dV. 
As can be seen in Eq. (18), this conductance only de- 
pends on the transverse components, i.e. the sum and 
difference of Xnm ^^'^ Xnmi since the term containing 
Xnm is proportional to nN -f m • M (= 0). Therefore, 
despite the presence of possible thermal noise, such a 
set-up would benefit from a very low current noise since 
most of the noise would be related to the spin fiuctua- 
tions, that is, the noise we want to measure. This can be 
seen by identifying the spin-dependent current operator 
with 5i{t) = TiS{t)-s, where s = ^p^,^, ^^crcka'- The 



current-current correlation function is then given bjP^'^ 
{Si{t)Si{t')) = T^s ■ (S(t)S(i')) • s, where we average over 
the dynamics of the localized spins and over the ensem- 
ble of the tunneling electrons. Under the condition that 
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nN + m • M = 0, the total dc current / is proportional 
to f_^s ■ (S(<)S(i')) ■ sd{t - t'), and since the shot 
noise is approximately (-^shot('^)) ~ the signal-to-noise 
ratio is about unity. 



V. CONCLUSIONS 

The theory we propose here for STM/STS mea- 
surements (essentially) con tain s the well-known TersofF- 
Hamann term,^°^^^^ Eq. (12), which maps the energy 
dependence electronic and magnetic densities of states, 
for energies defined by the Fermi level, of the substrate, 
and the applied bias voltage. For the situation where 
magnetic impurities are located on the substrate, our the- 
ory also contains a contribution which is proportional to 
the magnetic moment of the local spins, and allows for 
a quantitative analysis of spin moments using SP-STM. 
This is, hence, different from the currently most common 
analysis, which is solely focused on the energy depen- 
dent spin-polarization of the DOS. We show here that 
studies and analyses of both the dynamical conductance, 
dI{V,t)/dV or dI[V,u))/dV, and the total conductance 
dI{V)/dV can be linked to the dynamical, {S^{t)) or 
{S^{uj)), and total, J {S^{uj))dLd/2TT, magnetic moments, 
respectively, of the local spins. This contribution gen- 
erates different response at the local spins depending on 
their relative orientation compared to the spin moment 
on the SP-STM tip. To our knowledge, this contribution 
has not been discussed before, and the direct relation 
between the differential conductance and the local spin- 
moment is expected to make a significant impact on the 
potential capabilities in using SP-STM, and to extract 
quantitative information about local spin moments. 

We also point out that despite the present difficulties 
to experimentally record the time-dependent component 
in the conductance in a high GHz or even THz regime, 
it is important to observe that the dynamical component 
to the conductance opens new possibilities for the gen- 
eration of high frequency ac currents and/or voltages. 
Thus, by using local spin moments adsorbed onto metal- 



lic surface, such a scenario requires at least one mag- 
netic electrode which produces a local magnetic field that 
can interact with the local spin moment and, thus, pro- 
vide a time-dependent net contribution to the tunneling 
current, and conductance. Having in mind anisotropy 
fields from the substrate acting on the spin moment, 
it would be possible to generate spin-dependent ac cur- 
rents/voltages in the THz regime, since the anisotropy 
fields may be of the order of meV. 

It is important to note that our description goes be- 
yond the treatment reported in e.g. Refs. [23:24,. since 
we also include effects from the local spin-spin interac- 
tion between de-localized spin built up by the tunneling 
electrons, and the localized spin of the adatom. These 
spin-spin interactions are included already in the tunnel- 
ing matrix element, c.f. Eq. ([T]), and describe that the 
tunneling electrons of different spins are subject to dif- 
ferent tunneling barriers, i.e. spin-dependent tunneling 
barriers. 

The theory presented here for obtaining quantitative 
spin-resolved information in a SP-STM experiment pro- 
vides an alternative to the method used in Ref. 26, where 
a mean-field model based on the Weiss molecular fielcP^l 
was used. Future work, primarily of experimental nature, 
will judge which approach is the most reliable one. 
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